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Extended Abstract 

Problem description Traditional attitude control design of roosablc launch vcMclcs halves 
independent design of autopilot and control allocation modules [SBB99], [SK97] [SHJ^], JSO J 
[HodOOl Unfortunately, this results in the potential for overly aggressive commands m the autopilot 
resulting in a loss of performance due to actuator saturation, particularly if the autopilo may suffer 
from integrator wind-up [HH01]. This unfortunate situation can arise from, c.g., actuator pos^^^ 
limi ts that require that the actuator command vector J c (f) E A where A is the sc o 
actuator commands defined by 

A = {5 : 5- < 5 < (5+} 

where 5 and 5+ arc minimum and maximum position commands, respectively. Control allocation 
is limited by the attainable moment set T [BD95], [Dnr93], [Dur99], defined as 

T = {r : 35 E A and G6 = r} 

where G is the current Jacobian (control derivatives) matrix from the vehicle actuator condition 
vectort 6 to the vehicle body torques n 

„ a On 
G ~ls 

On-liuc calculation of the entire attainable moment set T is not a prMtical option for the following 
reasons: 

1. Actuator models (aerodynamic) arc approximate at best. 

2. Actuator failure will significantly modify the attainable moment set T. 

We propose instead to calculate a point-wise “snapshot” of T as shown in Figure 1. Given a torque 
command r c we compute a local attainable moment set T(f c ) dctcrmimnmg the maxmmm and 
minimum torque limits in each channel (x =roll, y =pitch, z =yaw) while holding the other torque 
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Figure 1: Attainable moment set snapshot 


values constant. The dimensions of this local attainable moment set can be calculated by a linear 
programming problem, c.g., given the current system Jacobian matrix G, 
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These torque limits may be of use in two scenarios: 

1. in communicating overall actuator torque limits to the autopilot and autocommandcr so that 
autopilot and/or guidance commands may be appropriately adjusted, and 

2. in flight scenarios where control allocation is required to divide torque commands a primary 
set of actuators (c.g, acrosurfaccs) and a secondary backup set of actuators (c.g., reaction 
thrusters). 


Torque limits computation and linear programming The computation of the local attain- 
able moment set can be posed as a set of six linear programming problems 

ma xJ(i) J{x) = c r x (0-1) 

subject to Ax = b 

x~ < x < x+ 
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where c T and A arc construct from appropriate rows of the Jacobian matrix G and the Rector b* 

the error Ax = b [Luc84].) Standard codes 1 available for the linear programming problem such 
dsplp and lp-solve arc inappropriate for this problem, since 

1. These codes arc design for large sparse problems, where as our A matrix is dense and w 
very few unknowns - at most 10 or 20. 

2. These algorithms can require several iterations to converge to an optimal solution, where as 
we require fast operation (low computational overhead). 

3 Our problem does not change too much from one time-step to the next (system . Jtmobim* 
do not change much over one sampling period), so our method should male use of earlier 

solutions. 

4. Further, our problem docs not require an exact optimal solution - just a good approximation. 

We therefore propose the use of a limited-iteration simplex method LP solution that uses resets 
(ac^e coSto nt set) of the previous iteration to compute initial values for the next iteratiom 
Kence indicates that two simplex iterations is sufficient for adequate ™ 

method requires the inversion of at most 24 2 x 2 matrices (2 for each limit, 2 limits PC r a ^h 

We show preliminary results of our method in Figure 2. The routine Ipsolve is a full simplex 
LP solver that exits upon convergence or detection of an infeasible problem (notice tj “ c s P‘ k 
t ~ 13scc) The routine lplter is an implementation of our fast LP solver. Both Ipsolve 
lplter have an outermost m-file script that is used to call a C-code implementation of om fast 
iteration. dmTrqLim is a full C-code translation of lplter that is included for comparison. Observe 
that, except when the input torque command is infeasible, lplter dose y tr s a ^ . 

torauc limits even though it is limited to at most two simplex iterations per time step. This ex P 

simulated actuators; similar results were obtained with much larger numbers 

of actuators in other simulation tests. ........ 

These torque limits can then be used to adjust autopilot control parameters on-hne ; so Ur* he 
anfrmilot can respond appropriately to either saturations or permanent failures. One potential 
method for using torque limits is presented explicitly in [SBB99]. We arc investigating t c use o 
this method in closed-loop control simulations of the X-33 launch vehicle. 
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